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ABSTRACT 


The  problem  of  extending  transport  property  measurements  to 
conditions  of  high  pressure  p  <  10^  atm  and  temperature  T  <  20,  000  K 
beyond  the  range  of  conventional  techniques  is  summarized.  The  ad¬ 
vantages  of  unconventional  techniques  over  conventional  transport 
property  measurement  techniques  in  regions  where  they  overlap  are 
presented  for  application  at  extreme  pressure-temperature  conditions. 
The  application  of  ultrasonic  absorption  and  sound  speed  to'  the  deter¬ 
mination  of  thermodynamic  and  transport  properties  at  high  pressures 
is  discussed.  Measurements  are  presented  at  300  K  and  pressures 
up  to  100  atm  in  the  argon- nitrogen  and  helium-argon  systems.  The 
schlieren  differential  interferometer  is  applied  to  thermal  conductivity 
measurements  at  high  pressures  and/or  temperatures.  Some  thermal 
conductivity  measurements  in  nitrogen  and  air  are  reported  at  tem¬ 
peratures  between  2000  and  5000  K. 
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SECTION  I 
INTRODUCTION 

Experiments  described  in  the  present  report  include  ultrasonic 
measurements  in  the  high  pressure  (p  <  400  atm)  region  as  well  as 
schlieren  interferometer  measurements  in  the  high  temperature 
(1000  <  T  <  5000  K)  region.  Ultrasonic  measurement  of  sound  ab¬ 
sorption  and  velocity  were  made  in  nitrogen,  argon,  two  mixtures 
of  nitrogen  and  argon,  and  a  50-50  helium-argon  mixture.  The 
pressure  of  these  measurements  ranged  up  to  400  atm  at  ~300  K. 
Schlieren  interferometric  measurements  of  thermal  conductivity  of 
nitrogen,  air  and  a  helium- argon  mixture  at  temperatures  up  to 
5000  K  and  pressures  about  1  atm  were  obtained.  The  experimental 
work  was  done  to  extend  the  information  on  transport  properties 
needed  for  engineering  design  at  high  pressures  and  temperatures  as 
well  as  for  long  range  support  of  theoretical  developments. 

Over  the  past  two  years  we  have  assessed  the  problem  areas 
in  transport  properties  of  gases'1'  and  designed  experimental  tech¬ 
niques  to  provide  data  in  these  areas.  The  problem  areas  we  have 
investigated  are  shown  graphically  in  Figure  1  for  air.  Also  shown 
are  techniques,  we  feel,  will  provide  useful  data  in  their  respective 
ranges.  Higher  order  collisions  at  high  densities,  diffusion  coef¬ 
ficients  under  conditions  where  reaction  conductivity  is  important, 
and  ionized  gases,  are  the  major  problem  areas  we  have  covered. 

In  addition,  at  low  pressures  (1  atm)  and  intermediate  temperatures 
(1000°C)  there  are  questions  (Refs.  3,4)  about  the  available  viscosity 
measurements  which  are  the  major  experimental  inputs  for  the  deter¬ 
mination  of  accurate  molecular  potential  functions.  The  potential 
functions  are  needed  in  engineering  practice  for  the  calculation  of 
thermodynamic  and  transport  properties. 

The  solid  lines  in  Figure  1  show  the  P-T  curve  along  which 
higher  order  collisions  cause  a  30%  and  70%  increase  in  thermal  con¬ 
ductivity  respectively  over  the  low  pressure  values.  Pressures  above 
and  temperatures  below  these  lines  result  in  density  effects  of  even 
greater  importance.  The  heat  transfer  measurement  provides  thermal 
conductivity  and  the  ultrasonic  measurements  provide  viscosity  and 
thermodynamics  (specific  heat  and  compressibility).  A  third  tech¬ 
nique  which  uses  schlieren  interferometric  determination  of  the  tem¬ 
perature  field  about  a  heated  or  cooled  cylinder  to  determine  thermal 

*We  have  not  considered  the  important  problems  which  arise  near  a 
phase  change  (Refs.  1,2). 
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conductivity  was  also  conceived  for  use  at  high  pressures  and  tem¬ 
peratures.  The  end  wall  heat  transfer  and  the  interferometric  method 
of  determining  thermal  conductivity  are  not  affected  as  much  as  the 
conventional  techniques  for  high  pressure  applications  by  convection 
and  radiation.  The  interferometric  measurement  is  capable  of 
operating  at  higher  temperatures  than  conventional  techniques.  In 
addition,  some  questions  have  arisen  about  the  validity  of  previous 
viscosity  measurements  which  do  not  apply  to  the  ultrasonic  technique. 
Finally,  thermodynamic  properties  (specific  heat,  compressibility)  of 
the  high  pressure  gas  which  are  not  always  known  accurately  at  the 
pressures  and  temperatures  of  interest  are  a  by-product  of  the  ultra¬ 
sonic  velocity  measurement. 

The  pressure-temperature  range  where  dissociation  is  im¬ 
portant  is  indicated  by  the  dotted  lines.  The  thermal  conductivity  of 
a  simple  dissociating  diatomic  gas  is  discussed  to  illustrate  the 
problems  involved.  The  equilibrium  thermal  conductivity  of  such 
a  dissociating  gas  is  given  by 


X  =  X 


+  pD  (h  -  h  ) 
m  a 


8  a 

a  t 


where  a  is  the  degree  of  dissociation,  hm  and  h&  are  the  enthalpy  of 
molecules  and  atoms,  respectively,  and  is  the  thermal  conduc¬ 
tivity  due  to  translational  energy  (frozen  thermal  conductivity), 

8  u 

p  D  (h  -  h  )  -r-=,  is  the  thermal  conductivity  due  to  diffusion  of  dis¬ 
til  a  d  1 

sociated  species  to  cooler  regions  with  subsequent  recombination 
(chemical  thermal  conductivity).  The  chemical  thermal  conductivity 
is  up  to  6  times  the  frozen  thermal  conductivity  depending  on  the  pres¬ 
sure  and  temperature  as  indicated  in  Figure  1.  Examination  of  the 
scanty  evidence  available  at  high  temperatures  from  the  Maecker  arc 
and  the  schlieren  interferometer  indicates  minimum  uncertainties 
(Ref.  7)  of  thermal  conductivity  of  50-100%  near  the  50%  dissociation 
temperature  in  nitrogen  and  air.  The  thermodynamic  quantities,  p  , 

8  c 

h&,  h^  and  may  be  calculated  accurately.  Accepting  an  uncer¬ 

tainty  in  the  diffusion  coefficient  (D)  of  20%  (as  is  often  quoted  in 
existing  tabulations)  the  large  enthalpy  associated  with  dissociation 
(hm-ha)  transforms  the.  20%  uncertainty  in  the  calculated  diffusion 
coefficient  directly  into  a  20%  uncertainty  in  thermal  conductivity. 
Accurate  determination  of  the  diffusion  coefficients  between  2000  and 
8000  K  is  an  immediate  engineering  need.  The  most  serious  need 
exists  around  4000  K  at  pressures  between  1/2-30  atm. 
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The  techniques  shown  in  Figure  1  for  application  in  the  dis¬ 
sociating  region  are  schlieren  interferometry  and  the  ultrasonic 
technique.  Schlieren  interferometry  has  been  used  during  the  past 
year  to  measure  equilibrium  thermal  conductivities  of  nitrogen, 
air  and  a  helium-argon  mixture  at  pressures  between  0.  1  and  1  atm 
and  at  temperatures  about  4000  K.  The  ultrasonic  technique  may  in 
some  cases  be  used  to  obtain  the  frozen  thermal  conductivity.  The 
chemical  thermal  conductivity  and  therefore  the  diffusion  coefficients 
may  in  principal  be  obtained  from  a  combination  of  the  schlieren  and 
ultrasonic  techniques. 

o 

The  high  temperature  range  above  10,  000  K  has  not  been 
investigated  under  the  present  contract  but  it  is  one  of  the  areas  of 
major  uncertainty.  The  transport  properties  of  ionized  gases  be¬ 
tween  0.  1  and  100  atm  are  an  engineering  requirement  for  the 
development  and  use  of  high  enthalpy  test  facilities.  The  techniques 
which  have  been  successfully  applied  to  this  area  are  the  Maecker  arc 
technique  (Refs.  5,  6,  7)  and  the  ultrasonic  technique  (Refs.  8,  9,  10). 
Electrical  conductivity  measurements  are  also  useful.  The  ultrasonic 
technique  requires  reaction  rates  for  ionization  and  collision  numbers 
for  electronic  excitation  to  obtain  transport  properties.  The  Maecker 
arc  technique  requires  electrical  conductivity  and  radiative  properties 
of  the  gas  to  be  known  before  the  thermal  conductivity  can  be  deter¬ 
mined.  The  schlieren  interferometric  technique  on  the  end  wall  of 

k  o 

the  shock  tube  may  be  used  to  obtain  q  between  8000  and  15,  000  K 

P 

with  no  knowledge  of  other  gas  properties  except  the  Gladstone- Dale 
constant  (Refs.  11,  12,  13). 

Summarizing,  the  methods  discussed  below  include  end  wall  heat 
transfer  and  ultrasonic  absorption  for  high  pressure  transport  property 
measurements.  These  experiments  are  complementary  since  end  wall 
heat  transfer  provides  thermal  conductivity  and  ultrasonics  provide 
viscosity  and  thermodynamics.  A  new  experiment  involving  schlieren 
measurements  of  the  density  field  about  a  horizontal  cylinder  in  con¬ 
vective  flow  is  mentioned  to  provide  high  density  thermal  conductivity 
at  high  temperatures  (T  ~  3000  K).  The  schlieren  interferometric 
measurement  in  the  end  wall  boundary  layer  is  also  presented  as  a 
means  of  measuring  equilibrium  thermal  conductivity  in  the  low  pres¬ 
sures,  high  temperatures  region  where  chemical  reactions  are  in 
progress.  This  measurement  combined  with  ultrasonic  measurements 
of  frozen  thermal  conductivity  provides  the  reaction  conductivity.  The 
reaction  conductivity  is  the  chief  source  of  pressure  dependence  in 
chemically  reacting  gases  (dissociating  air)  and  is  therefore  an  immediate 
engineering  requirement.  Finally,  new  measurements  relating  to  conduc¬ 
tivity,  viscosity,  diffusion  coefficients  and  thermodynamic  properties  of 
several  gases  are  presented. 
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SECTION  II 

THEORY  OF  ULTRASONIC  ABSORPTION  AT  HIGH  PRESSURES 


2.  1  CONSERVATION  EQUATIONS  AND  EQUATIONS  OF  STATE 


The  ultrasonic  method  for  determining  transport  properties  of 
gas  was  used  to  determine  the  effects  of  density  on  the  transport  and 
thermodynamic  properties  of  gases.  The  general  range  of  application 
of  the  ultrasonic  techniques  (including  Brillouin  scattering)  was  dis¬ 
cussed  in  last  year1  s  report  (Ref.  7). 


The  present  discussion  provides  working  equations  for  ultra¬ 
sonic  absorption  and  dispersion  which  are  valid  at  high  densities 
(p<  400  atm).  The  usual  expressions  for  sound  absorption  and  velocity 
in  low  pressure  gases  have  the  density  eliminated  in  favor  of  the  pres¬ 
sure  using  the  equation  of  state.  However,  density  is  the  fundamental 
variable  and  is  simpler  to  keep  at  high  densities.  The  potential  energy 
of  a  molecule  in  the  field  of  its  neighbors  becomes  an  important  contri¬ 
bution  to  the  internal  energy  at  high  densities.  This  leads  to  the  appear¬ 
ance  of  a  term  proportional  to  the  excess  density  in  the  equation  of  state 
for  the  internal  energy.  However,  the  interaction  of  a  molecule  with  its 
nearest  neighbor  is  not  sufficient  to  affect  the  rotational  and  vibrational 
(Ref.  14)  relaxation  times  up  to  densities  approaching  liquid  densities. 
The  change  in  relaxation  time  which  does  occur  is  due  to  a  change  in 
collisional  time  rather  than  the  transition  probability.  Finally,  an 
additional  effect  which  becomes  important  near  liquid  densities  is  a 
bulk  viscosity  due  to  higher  order  collisions. 


The  basic  phenomena  which  govern  the  sound  absorption  and 
speed  are  conservation  of 


energy 


momentum 


mass 


au 

+  £ 

a  v 

\ 

a  t 

=  0 

at 

p 

a  x 

P 

a  x 

8  v 

+  - 

1 

(t)’  + 

4  \  a2v 

—  r>  |  - 

at 

p 

a  x 

"  P 

3  7  a  2 
'  a  x 

rsj 

JLe_ 

at 

+  p 

0  V 

dx 

= 

0 

=  0 


(1) 

(2) 

(3) 
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and  the  equations  of  state  for  pressure 


fst 


JL 

p 

KT  p  -  P  T 

<V 

/  8U  \  ~ 

/dU\  ~ 

and  internal  energy  U  = 

(Tt)p  t  + 

(iSVt  p 

(4) 

(5) 


where  the  tilde  represents  the  fluctuating  part  of  the  dependent  variables 
and  the  dependent  variables  are  the  internal  energy  U,  temperature  T, 
density  p  ,  pressure  p,  and  velocity  v.  The  linearized  conservation 
equations  are  exact  for  typical  sound  waves  because  the  fluctuation  of 
temperature,  pressure,  and  density  in  a  typical  sound  wave  are  about 
10"'  of  the  ambient  values.  The  only  assumptions  involved  in  the  equa¬ 
tions  of  state  are  those  for  small  departures  from  equilibrium,  the 
isothermal  compressibility  K-p,  thermal  coefficient  of  expansion  p; 
specific  heat  Cv  =  (3U/8  T)p  ,  and  (8U/9p)T  are  essentially  the  equilib¬ 
rium  properties  of  the  gas. 

2.  2  SOUND  SPEED  AND  SPECIFIC  HEAT  RATIO  DETERMINATION 

The  assumption  that  the  ultrasonic  pulse  may  be  made  up  of  expo¬ 
nentially  damped  plane  waves  leads  to  an  expression  for  the  absorption 
coefficient  a  and  sound  speed  c.  The  viscosity  and  thermal  conductivity 
are  set  equal  to  zero  in  the  sound  speed  calculation  since  they  have  a 
negligible  effect  on  the  sound  speed  in  the  continuum  density  range 
(Refs.  9,  10).  The  condition  that  Eqs.  (1-5)  have  a  plane  wave  solution 
gives 


1_ 

2 


c 


C  + 
v 


P  Kr 


C 

v 


(6) 


Using  the  thermodynamic  identity 


(7) 
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the  sound  speed  becomes 


2 

c 


(8) 


The  sound  speed  measured  as  a  function  of  both  temperature  and 
pressure  is  sufficient  to  determine  all  the  thermodynamic  properties  of 
the  gas  (Ref.  7).  However,  in  the  present  instance  it  is  more  valuable 
to  use  compressibility  factors  from  pV  T  measurements  to  obtain  the 
specific  heat  ratio  y.  The  sound  speed  in  terms  of  the  compressibility 
is 

‘2  ■  [i-f  (|f)T]  1  <»> 

Equation  (9)  is  useful  (as  is  described  below)  to  obtain  specific  heat  ratios 
for  reducing  sound  absorption  to  thermal  conductivity. 

2.  3  THERMAL  CONDUCTIVITY  AND  VISCOSITY 

The  absorption  due  to  thermal  conductivity  and  viscosity  referred 
to  as  the  classical  absorption  ac  may  be  similarly  derived  from  Eqs.  (1-5) 
(Ref.  10).  The  resulting  expression  is 


<ji  c  L  p 


where  p  is  the  density,  c  sound  speed,  C^  specific  heat,  y  specific  heat 
ratio,  rj  viscosity  and  X.  thermal  conductivity.  The  viscosity  and  thermal 
conductivity  terms  increase  with  pressure.  However,  the  sound  speed 
cubed  in  the  denominator  of  (1)  increases  about  the  same  amount  between 
1-200  atm.  The  result  is  that  pec/u^  is  essentially  constant  in  this  pres¬ 
sure  range. 

The  thermodynamic  quantities  c,  y  ,  and  Cp  are  known  very  accu¬ 
rately.  The  thermal  conductivity  or  viscosity  may  therefore  be  determined 
when  independent  measurements  of  either  transport  property  are  available. 
On  the  other  hand,  the  measured  sound  absorption  may  be  used  to  test  the 
validity  theoretical  predictions  of  rj  and  y  with  no  other  measurements 
involved. 
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2.  4  BULK  VISCOSITY 

An  additional  transport  property,  the  translational  (or  intrinsic) 
bulk  viscosity  becomes  important  at  densities  approaching  the  critical 
density  p  c.  This  bulk  viscosity  should  not  be  confused  with  that  arising 
from  the  presence  of  internal  modes  which  is  also  present  in  polyatomic 
molecules.  The  translational  bulk  viscosity  may  be  considered  physically 
due  to  the  fact  that  correlations  between  the  position  of  molecules  die 
away  ir.  periods  longer  than  the  binary  collision  time.  There  is  there¬ 
fore  a  component  of  pressure  which  is  proportional  to  the  time  rate  of 
volume  change  and  which  lags  behind  the  pressure  associated  with  binary 
collisions. 

The  translational  bulk  viscosity  in  liquids  has  been  the  subject  of 
controversy.  However,  recently  the  experiment  (Refs.  15,  16,  17)  and 
theory  (Ref.  18)  have  been  shown  to  be  essentially  in  agreement.  In 
particular,  recent  studies  of  argon  (Ref.  16)  in  the  dense  gaseous  region 
have  shown  that  the  p  ^  variation  of  bulk  viscosity  predicted  the  Rice- 
Alanatt  theory  is  correct  with  the  experimental  error  -10%.  Following 
the  suggestion  of  Madigosky  the  bulk  viscosity  may  be  estimated  by  the 
Enskog  expression 


(ID 


where  the  Enskog  factor  x  is  taken  from  Rosenbaum  et  al.  (Ref.  19). 

The  bulk  viscosity  enters  the  sound  absorption  through  the  con¬ 
servation  of  momentum  Eq.  (2).  The  classical  absorption  is  modified 
by  replacing  4/3  rj  in  Eq.  (1)  by  4/3  rj  4-  r] 1  .  The  important  comparison 
for  the  present  discussion  is  between  the  excess  viscosity 


=  3  (l  ‘ 


(12) 


and  the  bulk  viscosity.  These  two  quantities  are  plotted  versus  pressure 
at  300°K  and  234.  6°K  in  Figure  2  for  argon.  For  temperatures  of  300  K 
and  above  the  Enskog  approximation  should  be  slightly  higher  (<  15%)  than 
the  true  value  of  the  bulk  viscosity.  The  bulk  viscosity  is  Jess  than  37% 
of  the  excess  shear  viscosity  up  to  pressures  of  about  100  atm.  The  bulk 
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viscosity  term  is  only  15%  of  the  excess  thermal  conductivity  term. 
Therefore,  other  measurements  of  viscosity  and  theoretical  estimates 
_  of  the  bulk  viscosity  may  be  used  to  obtain  the  excess  thermal  conduc¬ 
tivity  from  sound  absorption  measurements  up  to  about  200  atm  with  a 
maximum  uncertainty  of  2%.  A  better  approach  is  to  predict  the  excess 
viscosity,  thermal  conductivity  and  bulk  viscosity  theoretically  and  use 
the  ultrasonic  absorption  data  to  test  the  validity  of  the  theory.  Because 
of  the  bulk  viscosity,  the  ultrasonic  absorption  is  more  sensitive  to 
density  effects  than  other  transport  property  measurements. 

2.5  MIXTURES 

The  sound  loss  in  a  mixture  includes  the  usual  absorption  due 
to  the  thermal  conductivity  and  viscosity  of  the  mixture  plus  terms 
which  arise  from  diffusion.  Diffusive  separation  of  species  due  to 
the  temperature  and  pressure  gradients  in  the  sound  wave  leads  to 
an  additional  term  in  the  energy  equation.  The  theory  for  low  density 
(~1  atm)  gases  follows  vigorously  from  kinetic  theory  (Refs.  20,21). 
High  density  gases  should  be  treated  by  the  techniques  of  irreversible 
thermodynamics  (Ref.  22).  Neglecting  the  first  and  second  derivatives 
of  compressibility  with  respect  to  mole  fraction  the  diffusion  losses 
are  given  by  the  low  density  formula  (p  c^  replacing  p) 


P  a 


2  2 

2tt  f  x^x^7 


D 


(p’  +  ^ 


P  +  -2^-a 


■K° 


12 


(13) 


where  x  is  the  mole  fraction,  M  the  atomic  weight,  D12  the  binary  dif¬ 
fusion  coefficient  and  a  the  thermal  diffusion  ratio.  The  subscripts 
1  and  2  refer  to  the  two  species  (M2  <  Mj).  The  quantities  p  and  p1  are 
given  by 


(3 


m2-mi 

M 


(14) 


P' 


± .  1^1 
M  y  \7j-l 


(14a) 


The  7,  7  j,  and  72  are  the  specific  heat  ratios  of  the  mixture  and  the 
pure  gases  respectively. 
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The  term  involving  the  thermal  diffusion  factor  is  relatively  small 
compared  to  p  and  J3 '  in  cases  where  the  diffusion  losses  are  significant. 
Diffusion  losses  for  several  gas  mixtures  are  given  in  Table  1.  Mixtures 
of  hydrogen  and  helium  with  other  gases  generally  show  about  50%  of  the 
absorption  due  to  diffusion.  Sound  absorption  in  these  gases  may  be  used 
to  determine  diffusion  coefficients.  In  gases  such  as  air  and  argon-nitrogen 
mixtures  the  sound  absorption  due  to  diffusion  is  about  2%  of  the  total  ab¬ 
sorption.  The  sum  of  thermal  conductivity  and  viscosity  may  therefore  be 
determined  as  outlined  above  for  pure  gases.  The  2%  losses  may  be  cor¬ 
rected  using  theoretical  diffusion  coefficients  to  give  the  sum  of  viscosity 
and  thermal  conductivity  with  better  than  0.  1%  accuracy. 

The  density  behavior  of  the  transport  properties  of  mixtures  are 
not  known  very  well.  The  thermal  diffusion  term  may  be  calculated  from 
the  best  theoretical  estimate  since  its  contribution  is  small.  The  major 
factor  is  the  density  dependence  of  the  binary  diffusion  coefficient.  Binary 
diffusion  coefficient  measurements  at  high  densities  suggest  that  pDj 2  is 
more  sensitive  to  density  than  self- diffusion  coefficients  and  decreases  as 
density  increases  as  indicated  by  theory  (Refs.  25,26,27,28)  and  experi¬ 
ment  (Ref.  24). 


Table  1.  Contribution  of  various  loss  mechanisms  to  sound  absorption 
in  gas  mixtures  at  1  atm  with  1  MHz  sound  waves  (estimates) 


Mixture 

Temp. 

°K 

(ap/f\ 

x  10"6 

(ap/f2) 

V 

x  10"6 

(ap/f2)^, 
x  10"6 

(ap/f2^ 
x  10"6 

P 

P' 

^■aT 

y  I 

i_rQ  a  -m 

300 

0.  117 

0.  084 

0.  17  5 

-  1.  636 

-1.  636 

0.  153 

iie-Ar 

1000 

0.  134 

0.  109 

0.  228 

-1.  636 

He-N2 

350 

0. 0558 

0.  070 

0.  029 

0.  067 

-1.  500 

-1.  167 

0.  143 

Ar-N, 

300 

0.  056 

0.  106 

0.  042 

0.  009 

-0.  351 

-0.  685 

0.  024 

2 

1000 

0.  161 

0.  144 

0.  138 

0.  012 

-0.  351 

-0.  797 

0.  049 

H  -Ar 

300 

0. 0873 

0.  087 

0.  333 

-  1.  808 

-2.  141 

0.  096 

2 

1000 

0.  347 

0.  149 

0.412 

-1.  808 

-2.  178 

0.  121 

N_  -Cl 

350 

0.  093 

0.  105 

0.  093 

0.  001 

-0.  133 

-0.  133 

0.  087 

2  2 

1000 

0.  151 

0.  112 

0.  240 

0.  000 

-0.  133 

-0.  068 

0.  027 

N  -CCL 

300 

0.  367 

0.  105 

0.  183 

0.  002 

-  0.  444 

-0.  217 

0.  017 

2  2 

1000 

0.  688 

0.  166 

0.  760 

0.  000 

-  0.  444 

-0.  052 

0.  030 
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The  absorption  due  to  the  diffusion  coefficient,  therefore,  de¬ 
creases  with  increasing  pressure.  An  examination  of  Eq.  (13)  and  the 
sound  speed  versus  pressure  curves  shows  that  pa-p  decreases  about 
30%  in  100  atm  due  to  the  in  the  denominator.  There  may  be  an 
additional  decrease  of  about  10%  due  to  the  decrease  in  the  diffusion 
coefficient  in  the  same  pressure  range. 

2.  6  INTERNAL  MODES-RQTATIONAL  RELAXATION 

The  absorption  due  to  thermal  relaxation  at  high  densities  may 
be  derived  by  breaking  the  internal  energy  U  into  a  translational  part 
and  a  component  due  to  internal  modes.  Coupling  a  relaxation  equation 
(Refs.  10,  15)  to  the  conservation  equations  then  leads  to  the  correct 
expression  for  sound  absorption  due  to  internal  modes.  However,  it 
is  more  rigorous  at  high  densities  to  use  arguments  based  on  irreversible 
thermodynamics.  The  resulting  expressions  are  identical  to  the  usual 
low  density  expressions. 

The  only  question  remaining  is  the  effect  of  density  on  the  re¬ 
laxation  time.  The  main  event  leading  to  a  change  of  the  internal  state 
of  a  molecule  is  a  hard  binary  collision.  The  number  of  such  collisions 
per  second  should  be  directly  proportional  to  the  density.  The  product 
of  relaxation  time  and  the  density,  pT  ,  will  be  a  constant  at  all  densities 
if  the  probability  (transition  probability)  that  a  collision  will  result  in  a 
change  of  internal  state  is  constant.  This  assumption  is  plausible  from 
a  physical  point  of  view.  This  assumption  is  reasonable  because  the 
colliding  molecules  are  perturbed  by  a  hard  collision  much  more  than 
by  the  long  range  interaction  with  their  neighbors  up  to  high  densities. 
Measurements  of  vibrational  relaxation  times  (Refs.  14,  29)  show  that 
p  T  is  constant  at  densities  up  to  about  twice  the  critical  density.  The 
same  result  holds  up  to  at  least  100  psi  for  rotational  relaxation  in  N^ 
as  indicated  by  the  agreement  between  theory  and  experiment  in  Figure 
9.  The  absorption  due  to  rotational  relaxation  at  high  densities  reduces 
to 


where  r}Q  is  the  low  density  viscosity  and  R  is  the  gas  constant  per  gram. 
The  low  density  viscosity  enters  because  the  binary  collision  relaxation 
(Tc)  time  was  calculated  assuming  p  Tc  was  a  constant,  that  is 
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The  subscript  o  refers  to  the  low  density  limit.  The  density  sound  ab¬ 
sorption  product  indicated  by  Eq.  (15)  is  essentially  constant  between 
1-  100  atm. 

2.7  SUMMARY 

The  above  discussion  shows  how,  in  principle,  the  transport 
properties  of  gases  may  be  determined  at  high  pressures  using  ultra¬ 
sonics.  The  sum  of  thermal  conductivity  and  viscosity  may  be  deter¬ 
mined  from  Eq.  (10).  The  absorption  due  to  relaxation  of  internal 
modes  may  be  subtracted  from  the  total  absorption  to  obtain  the  sum 
of  thermal  conductivity  and  viscosity  in  diatomic  and  triatomic  gases. 
The  frequency  may  be  chosen  sufficiently  high  so  that  all  internal 
modes  except  rotation  contribute  less  than  0.  1%.  Equation  (15)  then 
is  typically  30%  correction  term.  Finally,  ultrasonic  absorption  in 
mixtures  may  be  used  to  measure  diffusion  coefficients  or  the  sum  of 
viscosity  and  thermal  conductivity  depending  on  the  difference  in 
molecular  weights  of  the  species  involved.  Diffusion  losses  in  mix¬ 
tures  of  helium  and  hydrogen  with  heavier  gases  are  ~50%  of  the  total. 
Therefore,  absorption  and  diffusion  coefficients  of  such  mixtures  may 
be  obtained.  On  the  other  hand,  diffusion  losses  in  mixtures  of  other 
gases  are  small  (typically  less  than  ~1%)  and  the  sum  of  viscosity  and 
thermal  conductivity  may  be  determined  as  described  above. 
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SECTION  in 

APPLICATION  OF  THE  ULTRASONIC  TECHNIQUE  TO 
HIGH  PRESSURE  GASES 

3.  1  PULSE  ECHO  THROUGH  TRANSMISSION  TECHNIQUE 

High  pressure  cells  suitable  for  measuring  ultrasonic  absorption 
and  dispersion  have  used  both  movable  probes  (Ref.  30)  and  fixed  probes 
(Ref.  31).  Movable  probe  apparatus  yields  curves  of  amplitude. versus 
position  and  transit  time  versus  distance  directly.  The  fixed  probe  tech¬ 
nique  must  use  echoes  to  obtain  more  than  one  path  length.  In  addition, 
corrections  for  the  diffraction  losses  and  transmissivity  of  the  ultrasonic 
mirrors  must  be  made.  The  advantages  of  the  fixed  probe  technique  are 
due  to  its  mechanical  simplicity.  There  are  no  errors  associated  with 
variation  of  alignment  with  position.  There  are  no  rotating  shafts  to  seal 
against  high  pressure.  Finally,  the  mechanical  simplicity  of  the  fixed 
path  technique  is  very  important  at  high  temperatures  (500  C)  and  pres¬ 
sures  (8000  psi). 

The  fixed  path  probes  used  in  the  present  experiment  are  shown 
in  Figure  3.  The  details  of  the  pressure  cell  are  not  shown.  The  probes 
are,  however,  attached  to  the  cover  and  the  floor  of  a  steel  pressure 
bomb  with  suitable  high  pressure  inlets  for  gas,  electricity,  internal 
water  cooling,  thermocouples,  and  coaxial  lines  for  the  rf  pulses  which 
drive  the  ultrasonic  probes.  The  pressure  cell  used  with  the  probes  is 
internally  cooled  so  that  temperatures  up  to  ~2000°K  at  pressure  up  to 
200  atm  could  be  used.  This  cell  (referred  to  as  I)  was  used  to  obtain 
the  sound  absorption  measurements  reported  below.  A  second  cell 
(referred  to  as  II)  externally  heated  (built  by  Aminco)  which  will  go  to 
1000  atm  at  300°K  and  400  atm  at  800°K  was  used  to  obtain  the  sound 
speed  reported  between  100  atm  and  400  atm.  The  probes  in  cell  II 
were  similar  to  those  shown  in  Figure  3  except  that  both  probes  were 
mounted  on  the  cover  of  the  pressure  cell. 

The  major  problems  associated  with  acoustic  measurements 
in  gases  at  high  temperature  and  pressure  are  alignment,  acoustic 
isolation,  rf  transmission,  ringing  and  diffraction  corrections.  The 
probes  are  held  in  alignment  by  the  stainless  steel  tubing  which  also 
serves  as  the  ground  for  the  rf  transmission  lines.  Ceramic  insulators 
allow  the  coaxial  line  to  be  used  at  high  temperatures.  The  probes  are 
roughly  aligned  by  maintaining  the  mechanical  parts  as  true  as  possible. 
However,  final  alignment  must  be  adjusted  while  examining  the  first 
pulse  and  the  echo  by  an  oscilloscope.  The  previous  systems  mentioned 
above  used  a  solid  tube  for  alignment. 


12 


AEDC.TR.69.78 


The  sound  energy  is  more  easily  transmitted  through  the  trans¬ 
ducer  mount  into  the  walls  of  the  pressure  cell  than  through  the  desired 
gas  path.  This  is  because  of  the  large  mismatch  between  the  gases  and 
the  solids  in  the  low  pressure  range  (p  <  100  atm).  Improper  acoustic 
isolation  leads  to  sound  transmission  through  the  walls  which  overshadows 
the  desired  signal  through  the  gas.  The  only  technique  which  provides 
sufficient  isolation  between  the  cell  walls  and  the  transducer  under  the 
conditions  of  the  present  experiment  is  minimum  mechanical  contact  be¬ 
tween  the  receiver  and  the  transmitter.  Mounting  the  probe  holder  on 
the  coaxial  line  as  shown  in  Figure  3  provides  minimum  mechanical 
contact. 


The  remaining  features  of  probe  design  are  best  understood  in 
terms  of  the  through  transmission  multiple  echo  technique  employed  in 
the  present  experiments.  The  photographs  in  Figure  4  show  a  series 
of  ultrasonic  pulses.  The  first  pulse  in  each  photograph  is  the  trans¬ 
mitted  pulse  as  it  encounters  the  receive  transducer.  The  following 
pulses  arise  by  multiple  reflections  in  the  gas  cavity  between  the  trans¬ 
mitter  probe  and  the  receive  probe.  Thus,  the  distance  between  the 
first  pulse  and  the  first  echo  is  twice  the  probe  separations.  There  is 
a  "tail"  on  each  pulse  which  interferes  with  the  following  pulse  at  close 
separations.  The  tail,  referred  to  as  ringing,  is  due  to  acoustic  energy 
trapped  in  the  transducers.  The  ringing  may  be  damped  for  a  crystal 
radiating  into  air  by  backing  the  transducer  with  a  lossy  material.  The 
damping  material  used  in  the  present  experiment  was  a  short  circuited 
PZT-5  ceramic  transducer. 

The  attenuation  of  the  ultrasonic  wave  is  given  by  the  natural 
logarithm  of  amplitude  ratio  of  the  first  pulse  and  the  first  echo, 
divided  by  the  path  length.  The  measured  attenuation  is  due  to  ab¬ 
sorption  during  reflection,  diffraction  losses  and  absorption  due  to 
transport  and  relaxation  phenomena.  The  attenuation  coefficient  due 
to  transmission  at  reflecting  surfaces  may  be  neglected  at  pressures 
up  to  100  atm  because  the  gas  acoustic  impedance  is  small  compared 
to  the  solid.  The  effective  attenuation  coefficient  due  to  reflection 
losses  is: 


where  p  cs  are  the  products  of  the  sound  speed  and  the  density  of  the 
gas  and  the  solid,  respectively.  The  absorption  coefficient  due  to  re¬ 
flection  losses  in  the  present  experiment  which  ranges  between  5  x  10“^ 
to  5  x  10-^  nepers/cm  is  negligible. 
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The  diffraction  losses  depend  only  on  the  wavelength  of  the  ultra¬ 
sonic  wave,  the  diameter  of  the  radiator  and  the  path  length.  The  sim¬ 
plest  loss  of  this  type  is  due  to  beam  spreading.  That  is,  the  energy 
diverges  so  that  a  certain  fraction  of  the  transmitted  beam  does  not 
intercept  the  receiver.  Beam  spreading  losses  amount  to  about  0.  01 
nepers/cm  under  the  worst  conditions  of  the  present  experiment. 

The  term  "diffraction  losses"  (Ref.  32)  is  more  commonly  used 
to  describe  apparent  losses  due  to  deviations  from  constant  pressure 
across  the  beam  which  are  inherent  in  the  piston  radiator  field.  In  the 
MHz  frequency  range  for  typical  radiator  dimensions  (a  ~1  cm)  diffrac¬ 
tion  losses  are  small  compared  to  the  absorption  due  to  transport 
phenomena.  However,  due  to  the  mounting  of  the  transducer  in  a 
slight  depression,  diffraction  losses  are  much  larger  than  usual. 
Typically,  a^ff  is  about  0.  2  nep/cm  in  the  present  experiments.  A 
flush  mounted  transducer  has  lower  diffraction  losses  but  is  more  dif¬ 
ficult  to  electrically  ground. 

3.2  DATA  REDUCTION 

The  separation  of  the  attenuation  due  to  diffraction  from  the 
absorption  due  to  transport  and  relaxation  effects  may  be  accomplished 
many  ways.  The  difference  between  the  measured  low  pressure  (~1 
atm)  attenuation  in  the  pressure  cell  and  the  theoretical  absorption 
could  be  used  to  obtain  the  diffraction  losses.  The  method  adopted  in 
the  present  experiment  takes  advantage  of  the  pressure  dependence  of 
the  various  loss  mechanisms.  As  mentioned  in  the  theory  section,  the 
classical  absorption  varies  as  l/p  .  The  absorption  due  to  diffraction 
losses  and  reflection  losses  is  essentially  independent  of  density.  There¬ 
fore,  a  plot  of  attenuation  versus  l/p  should  be  a  straight  line 

a  =  —  +  B  (18) 

P 

The  constant  B  is  the  diffraction  loss.  A  check  on  the  validity  of  this 
procedure  is  provided  by  the  fact  that  A  should  be  equal  to  the  density 
times  the  classical  absorption  at  atmospheric  pressure.  The  linear 
relation  clearly  holds  as  may  be  seen  from  Figure  5.  In  addition,  the 
coefficient  of  the  l/p  term  is  equal  to  the  classical  absorption  calculated 
from  measured  values  of  viscosity  and  thermal  conductivity. 

It  should  be  noted  that  this  technique  applies  over  a  limited  pres¬ 
sure  range.  The  assumptions  are  that  reflection  losses  (a  )  and  the 
sound  speed  (or  wavelength)  do  not  change  much  over  the  pressure  range 
of  the  data. 
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3.  3  SUM  OF  THERMAL  CONDUCTIVITY  AND  VISCOSITY 
MEASUREMENTS  IN  THE  ARGON-NITROGEN  SYSTEMS 

Measurements  of  sound  absorption  and  velocity  in  argon,  nitrogen 
and  argon- nitrogen  mixtures  at  pressures  up  to  100  atm  illustrate  the  type 
of  information  which  has  been  obtained.  The  P&c  versus  pressure  is 
shown  in  Figure  6  for  argon.  The  specific  heat  ratio  and  sound  speed 
from  the  present  measurements  were  used  with  independent  thermal 
conductivity  and  viscosity  and  Cp  measurements  to  calculate  p  a  c>  The 
agreement  between  the  various  experiments  is  excellent  over  the  whole 
pressure  range.  The  pressure  effect  on  viscosity  and  thermal  conduc¬ 
tivity  over  the  pressure  range  of  100  atm  amounts  to  about  25%.  The 
bulk  viscosity  contributes  about  10%  of  the  excess  or  2%  of  the  total  ab¬ 
sorption  in  the  gases  measured. 

The  normalized  absorption  (p  a  )  measured  in  nitrogen  is  shown 
in  Figure  7.  Nitrogen  has  the  additional  feature  of  rotational  relaxation. 
The  classical  absorption  is  shown  by  the  solid  line  in  Figure  7.  The 
absorption  due  to  rotational  relaxation  was  calculated  using  a  rotational 
collision  number  Zrot  =  4.  25,  determined  from  ultrasonic  measurements 
at  atmospheric  pressure.  The  total  absorption  passes  through  the  center 
of  the  data  over  the  entire  pressure  range. 

The  sound  absorption  versus  pressure  in  the  argon- nitrogen 
mixtures  was  difficult  to  handle  because  of  the  uncertainties  associated 
with  the  properties  of  mixtures  at  high  pressures.  The  transport  prop¬ 
erties  of  mixtures  may  have  a  different  pressure  behavior  than  pure 
gases  because  of  the  pressure  dependence  of  the  diffusion  coefficients. 
However,  the  technique  of  fitting  the  normalized  absorption  to  a  linear 
function  of  pressure  was  applied.  The  procedure  is  justified  by  the  fact 
that  the  coefficient  of  l/p  does  come  out  to  be  equal  to  the  known  absorp¬ 
tion  at  atmospheric  pressure.  The  rotational  collision  number  was  de¬ 
termined  for  the  mixture  from  atmospheric  pressure  measurements  and 
Eq.  (15).  The  sum  of  viscosity  and  thermal  conductivity  in  terms  of 
p  a  c  was  added  to  the  rotational  contribution  to  compare  with  the  mea¬ 
sured  p  a  .  The  results  are  displayed  in  Figures  8  and  9  along  with  the 
theoretical  total  absorption.  It  should  be  pointed  out  that  the  sum  of 
viscosity  and  thermal  conductivity  would  be  obtained  by  subtracting  the 
rotational  losses  from  the  total  absorption.  Accurate  specific  heat  data, 

7  and  C  ,  necessary  to  reduce  p  a  c  to  transport  properties  are  not  avail¬ 
able.  Inis  shortcoming  will  be  removed  by  measurements  of  sound  speed 
versus  pressure  at  several  temperatures  above  300  K. 


15 


AEDC-TR-69-78 


3.  4  SUM  OF  THERMAL  CONDUCTIVITY  AND  VISCOSITY 
MEASUREMENTS  IN  THE  HELIUM-ARGON  SYSTEMS 

Absorption  versus  pressure  in  a  0.  5  helium-0.  5  argon  mixture 
is  shown  in  Figure  10.  These  measurements  were  suggested  by 
J.  Sengers  for  comparison  with  some  high  density  results  on  diffusion 
coefficients  in  a  Lorentz  gas  (Refs.  26,  27,  28).  The  absorption  due 
to  diffusion  is  about  one-half  of  the  total  absorption.  The  classical 
absorption  shown  in  Figure  10  was  constructed  by  extrapolating  p  a  c 
from  low  pressure  values.  The  contribution  of  diffusion  was  calculated 
by  extrapolating  one  atmosphere  diffusion  coefficients  theoretically 
(Ref.  28)  to  high  pressures.  There  are  some  indications  (Ref.  33) 
that  the  thermal  diffusion  factor  could  increase  by  a  factor  of  3  at  high 
pressures.  This  effect  was  included  in  the  calculation  of  diffusion 
losses.  The  total  absorption  coefficient  thus  calculated  agrees  re¬ 
markably  well  with  the  measured  values. 

The  agreement  between  theory  and  experiment  in  Figure  10 
shows  that  the  theory  is  fundamentally  sound.  However,  the  usefulness 
of  the  ultrasonic  technique  for  the  determination  of  diffusion  coefficients 
at  high  pressures  is  limited  by  several  factors.  First,  the  above  cal¬ 
culations  show  that  the  size  of  the  diffusion  losses  decreases  rapidly  as 
pressure  increases  due  to  the  c^  in  the  denominator  of  Eq.  (13).  Second, 
the  thermal  diffusion  coefficient  becomes  more  important  at  high  pres¬ 
sures.  For  He-CO  mixtures  the  *Y  ~  1  a  factor  increases  to  about 

c.  y  T 

60%  of  the  diffusion  term  in  Eq.  (13)  at  100  atm. 

3.  5  SPECIFIC  HEAT  RATIO  DETERMINATIONS  FROM 
SOUND  SPEED  MEASUREMENTS 

Sound  speed  measurements  were  obtained  simultaneously  with  the 
absorption  measurements.  The  temperature  was  measured  with  two 
thermocouples  (see  Figure  8).  The  sound  speed  in  this  case  may  be 
used  to  determine  the  thermodynamic  properties  of  the  gas.  Specific 
heat  ratios  were  obtained  from  compressibility  measurements  (Ref.  34) 
and  the  sound  speeds  for  argon  and  nitrogen.  Specific  heat  ratios  for 
the  argon-nitrogen  mixtures  were  derived  based  on  the  assumption  that 
the  excess  viral  coefficients  were  small  (Ref.  35).  Specific  heat  ratios 
were  not  determined  for  the  helium-argon  mixtures  because  the  excess 
viral  coefficients  may  be  large  (Ref.  3  5). 

The  major  part  of  the  pressure  dependence  of  the  sound  speed 
comes  from  the  specific  heat  ratio.  The  factor 
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in  Eq.  (9)  differs  from  1.  00  at  most  by  2%  in  the  1-100  atm  range  for 
the  argon-nitrogen  system.  Thus  very  accurate  specific  heat  ratios 
may  be  obtained  from  rough  compressibility  measurements  and  ac¬ 
curate  sound  speed  measurements. 

A  schematic  diagram  of  the  electronics  used  for  the  sound  speed 
measurements  is  shown  in  Figure  11.  The  first  pulse  and  the  first  echo 
are  superimposed  on  a  dual  beam  oscilloscope  using  delayed  triggers. 

The  delayed  scope  trigger  is  obtained  by  putting  the  transmitter  trigger 
through  a  two-channel  digital  delay  generator.  The  oscilloscope  traces 
are  delayed  to  cancel  the  transit  time  of  the  respective  pulses  to  within 
0.  1  |J.sec.  Since  the  transit  times  are  typically  150  p.  sec  this  procedure 
gives  a  precision  of  better  than  0.  1%.  Delays  inherent  in  the  triggers 
and  errors  in  the  ultrasonic  path  length  were  corrected  with  calibration 
in  argon  and  helium  at  1  atm.  A  typical  oscilloscope  trace  is  shown  in 
Figure  12. 

The  sound  speed  versus  pressure  curves  for  the  argon-nitrogen 
system  is  shown  in  Figure  13.  Where  data  were  available  the  present 
data  are  in  agreement  with  other  work  to  within  0.  5%.  The  tempera¬ 
ture,  pressure  and  sound  speed  are  tabulated  in  Tables  2  and  3.  Also 
shown  are  values  of  specific  heat  ratio  7  derived  from  the  sound  speed 
and  compressibility  measurements  using  Eq.  (9).  The  mixtures  again 
were  a  problem  because  compressibilities  were  not  available.  Further 
measurements  at  different  temperatures  are  in  progress  to  obtain 
thermodynamic  properties  from  the  sound  speed  measurements  alone. 

3.  6  SUMMARY 

A  cell  for  measuring  sound  absorption  and  sound  speed  at  pres¬ 
sures  up  to  1000  atm  and  temperatures  up  to  800  K  has  been  built  and 
tested.  Attenuation  measurements  have  been  made  in  this  apparatus  with 
better  than  1%  precision.  There  are  large  diffraction  losses  (~0.  2  nep/cm) 
inherent  to  the  acoustic  probes  used  due  to  the  recessed  mounting  of  the 
acoustic  transducer.  Since  the  diffraction  losses  had  to  be  subtracted 
from  the  total  attenuation  coefficient,  the  precision  of  the  absorption  mea¬ 
surements  reported  was  reduced  as  much  as  15%  at  the  highest  pressures. 
While  these  experiments  were  in  progress  a  new  flush  mounting  technique 
was  developed  reducing  the  diffraction  losses  by  a  factor  of  10  (to  0.  02 
nep/cm).  This  increase  in  precision  will  be  reflected  in  future  absorp¬ 
tion  measurements. 
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Measurements  of  absorption  in  the  argon-nitrogen  system  show 
that  the  experiment  and  theory  are  in  agreement.  Thermal  conductivities 
could  be  obtained  from  the  smoothed  sound  absorption  data  and  other  vis¬ 
cosity  measurements.  However,  more  accurate  measurements  are  in 
progress.  The  measurements  of  sound  absorption  in  the  helium-argon 
mixture  show  that  diffusion  coefficients  are  difficult  to  obtain  from  ab¬ 
sorption  measurements  at  high  pressures.  This  is  due  to  the  decrease 
of  the  diffusion  losses  with  an  increase  in  pressure.  On  the  other  hand, 
the  decrease  of  significant  diffusion  losses  means  that  accurate  measure¬ 
ments  of  thermal  conductivity  and  viscosity  may  be  obtained  from  ab¬ 
sorption  measurements  in  gas  mixtures  at  high  pressures. 

Finally,  sound  speed  versus  pressure  measurements  were  made 
in  the  argon- nitrogen  and  argon-helium  systems.  The  major  reason 
for  sound  speed  variation  with  pressure  below  100  atm  is  the  change 
in  the  specific  heat  ratio;  therefore,  specific  heat  ratios  were  derived 
from  the  sound  speed  measurements  in  the  argon-nitrogen  system. 
Because  the  unknown  excess  viral  coefficients  are  large  in  the  helium- 
argon  mixtures  no  attempt  was  made  to  obtain  specific  heat  ratios  for 
these  gases.  Variation  of  pressure  above  100  atm  and  variation  of 
temperature  will  provide  both  compressibility  and  specific  heat  ratio 
from  sound  speed  alone.  Thermodynamic  properties  of  gas  mixtures 
at  pressures  above  200  atm  are  not  well-known;  therefore,  the  sound 
speed  measurement  is  very  useful. 
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Table  2.  Sound  speed  data  in  the  argon-nitrogen  systems 


Nitrogen 


0-200 

psi  Cell  I 

0- 

30,  000  psi 

Cell  II 

Sound 

Sound 

Pressure 

Temp. 

Speed 

Pressure 

Temp. 

Speed 

P 

T 

C 

7 

P 

T 

C 

7 

atm 

°K 

m/sec 

atm 

°K  • 

m/sec 

3.  93 

300.  3 

353.  0 

1.401 

8.  84 

302.  6 

349.  2 

1.  361 

3.  93 

299.  8 

353.  2 

1.405 

20.  41 

302.  5 

357.  6 

1.42  5 

14.  67 

299.  6 

355.  2 

1.422 

40.  82 

302.  6 

360.  9 

1.  449 

21.  54 

299.  6 

356.  3 

1.428 

61.  22 

3  02.  7 

366.  4 

1.  493 

28.  35 

299.  8 

358.  3 

1.  448 

81.  29 

303.  3 

371.  6 

1.  505 

34.  88 

299.4 

359.  7 

1.455 

102.  04 

303.  5 

378.  1 

1.  533 

41.  61 

299.  0 

360.  3 

1.463 

122. 44 

303.  5 

383.  9 

1.  556 

55.  22 

298.4 

362.  9 

1.484 

142. 86 

303.  2 

395.  4 

62.  36 

298.  6 

365.  2 

1.  504 

163. 37 

302.  7 

407.  6 

68.  62 

298.  6 

367.  3 

1.  509 

184.  0 

302.  3 

411.  9 

82.  90 

299.  0 

372.  5 

1.  534 

225.  5 

301.  3 

445.  9 

88.  76 

297.  8 

374.  6 

1.  550 

264.  3 

301.  9 

464.  2 

88.  62 

297.  8 

372.  2 

1.  531 

295.  2 

301.  3 

488.  0 

95.  69 

297.  8 

375.  9 

1.  553 

328.  6 

301.  3 

505.  7 

102. 36 

297.  8 

378.  4 

1.  565 

3  57.  8 

301.  2 

527.  0 

102.  2-2 

297.  8 

378.4 

1.  565 

50  N2-50  Ar 

P 

T 

C 

7 

P 

T 

C 

7 

7.  75 

295.  6 

327.  1 

68.  03 

306.  1 

339.  9 

14.  69 

296.  5 

328.  6 

81.  63 

305.  9 

342.  6 

21.  43 

299.2 

330.  2 

95.  24 

306.  0 

349.  9 

28.  23 

297.  7 

331.  7 

109.  5 

305.  6 

355.  5 

35.  44 

297.7 

333.  1 

122.  8 

305.  0 

360.  0 

41.  70 

297.  8 

334.4 

49.  18 

297.  7 

334.  9 

55.  24 

297.  7 

339.  2 

61.  43 

297.7 

338.  7 

68.  37 

297.  7 

339.  2 

74.  97 

297.7 

342.  0 

80.  00 

298.  2 

343.  0 
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Table  2.  Sound  speed  data  in  the  argon- nitrogen  systems  (contr  d) 


0-200 

psi  Cell  I 

Argon 

7 

0- 

30,  000  psi 

Cell  II 

7 

Pressure 

P 

atm 

Temp. 

T 

°K 

Sound 

Speed 

C 

m/sec 

Pressure 

P 

atm 

Sound 
Temp.  Speed 

T  C 

K  m/sec 

7.  35 

299.  8 

320.  6 

1.  648 

60.  88 

304.  9 

328.  2 

1.  787 

21.  91 

297.  6 

322.  0 

1.  686 

81.29 

304.  8 

331.4 

1.  840 

28.30 

297.  7 

322.9 

1.  708 

101.4 

304.4 

340.  2 

1.  961 

34.  22 

297.  7 

323.  8 

1.720 

121.  8 

304.4 

343.  6 

2.  018 

42.  25 

297.7 

324.  8 

1.  779 

140.  5 

303.2 

350.  5 

48.  64 

297.  7 

326.  9 

1.  811 

163.  3 

303.  2 

356.  5 

54.  83 

297.  7 

326.  6 

1.  808 

181.  6 

303.  2 

368.  1 

61.  84 

298.  0 

328.  1 

1.  844 

204.  8 

303.  9 

377.  2 

68.  10 

298.  0 

329.  5 

1.  849 

224.  8 

303.  9 

384.  5 

74.  70 

297.  8 

330.2 

1.  863 

245.  9 

304.2 

400.  4 

81.  78 

298.  2 

332.  2 

1.  890 

266.  3 

304.  6 

407.  3 

88.44 

298.  2 

333.  9 

1.  915 

289.  1 

305.  6 

418.  8 

95.24 

298.  2 

335.  6 

1.  940 

311.  6 

305.  4 

433.  3 

102. 04 

298.  3 

337.4 

1.  968 

329.  9 

305.4 

443.  4 

345.  6 

304.  4 

453.  5 

0. 

9  Ar-0. 

1N2 

P 

T 

C 

T 

7.  96 

300.4 

323.  8 

14.35 

301.  0 

323.8 

31.  64 

301.  0 

325.  0 

28.  57 

301.  0 

326.  1 

35.  24 

301.  0 

327.  6 

42.  38 

301.  0 

328.  5 

48.  98 

301.  0 

329.  0 

55.  24 

301.  0 

330.  4 

61.  50 

301.  0 

331.2 

68.  64 

301.  0 

332.  1 

75.  10 

301.  2 

334.  1 

81.  50 

301.  0 

334.  8 

82.  30 

301.  0 

336.  6 

95.  10 

301.  2 

338.4 

101.  91 

301.  2 

340.4 
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Table  3.  Sound  speed  data  in  a  0.474  argon- 0.  526  helium  mixture 


0-200 

psi  Cell  I 

Cell  II 

Pressure 

Sound 
Temp.  Speed 

Pressure  Temp. 

Sound 

Speed 

P 

atm 

T  C 

°K  m/sec 

7 

P  T 

r  0 

atm  K 

C 

m/sec 

7 

7.  35 

299.  0 

426.4 

14.  69 

299.  0 

428.4 

21.  77 

299.  0 

430.4 

28.  16 

296.  8 

430.  6 

35.  10 

297.4 

432.  1 

41.  77 

297.4 

434.  8 

41.  77 

296.4 

432.  5 

48.  44 

297.4 

436.  1 

49.  52 

296.  6 

436.  7 

54.  97 

297.  7 

430.  8 

61.  77 

297.  7 

441.  0 

68.  98 

297.7 

443.  1 

75.  10 

297.  9 

445.3 

81.  90 

297.  9 

447.  1 

89.  12 

297.  9 

447.  5 

96.  26 

297.  9 

449.  3 

102.  04 

297.  9 

451.  1 
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SECTION  IV 

SCHLIEREN  INTERFEROMETER 

4.  1  WOLLASTON  PRISM  DIFFERENTIAL  INTERFEROMETER 

A  differential  interferometer  measures  the  index  of  refraction 
gradient.  This  in  turn  is  directly  proportional  to  the  density  gradient 
for  many  systems  of  interest.  The  measurement  of  thermal  conduc¬ 
tivity  by  conventional  techniques  requires  temperature  derivatives  to 
be  calculated  from  temperature  measurements.  Since  the  differential 
interferometer  measures  gradients  directly  in  the  gas  it  has  some  in¬ 
teresting  applications  for  thermal  conductivity  measurements.  The 
problem  is  to  find  flows  where  the  temperature  gradients  in  the  con¬ 
servation  of  energy  equation  may  be  eliminated  in  favor  of  the  density 
gradient.  The  temperature  field  about  a  body  in  free  convection  and 
the  end  wall  boundary  layer  arjp  two  such  applications. 

The  "schlieren"  or  Wollaston  prism  differential  interferometer 
is  a  very  useful  device  in  its  own  right.  It  is  simple  and  not  as  sensitive 
to  vibration  as  the  usual  interferometer.  Wollaston  prism  differential 
interferometers  have  been  described  by  Lamb  and  Schreiber  (Ref.  6) 
for  arc  applications  and  by  Oertel  (Ref.  36)  for  shock  tunnel  applications. 
The  primary  interest  in  this  instrument  is  based  on  Smeets  (Refs.  37,  7) 
thermal  conductivity  determinations  using  the  end  wall  boundary  layer. 

A  schematic  of  the  basic  apparatus  is  shown  in  Figure  14.  The 
light  from  any  point  in  the  FX12  mercury  light  source  is  polarized  so  the 
electric  vector  is  45°  to  the  optical  axis  of  the  Wollaston  prism.  The 
prism  is  constructed  of  two  optically  active  quartz  wedges  (Ref.  6) 
cemented  together  with  crossed  optical  axis.  The  Wollaston  prism 
breaks  the  polarized  light  into  two  rays  of  equal  intensity  diverging 
at  an  angle  e.  The  lens  Lj  is  focused  through  the  Wollaston  prism 
and  the  polarizer  into  the  center  of  the  light  source.  Therefore  the 
diverging  rays  from  the  prism  are  converted  into  parallel  rays  by  Lj. 

In  addition,  the  rays  are  separated  by  a  distance  d.  A  similar  system 
recombines  the  two  rays  on  the  film  plane  of  the  camera.  The  lens 
is  focused  on  an  imaginary  object  in  the  center  of  the  test  region. 

When  the  test  region  is  uniform  a  fringe  pattern  develops  due  to 
phase  shifts  within  the  two  Wollaston  prisms  (Ref.  6).  The  location  of 
the  fringes  (or  equivalently  the  phase  shifts  in  the  system)  depends  on 
the  wavelength  of  the  light.  Thus  red  fringes  do  not  form  at  the  same 
place  on  the  photographic  plane  as  the  yellow  or  blue  fringes.  Visually 
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the  interference  pattern  from  white  light  separates  the  colors.  On  a 
photograph  the  effect  is  to  show  fuzzy  fringes  since  the  camera  does 
not  discriminate  among  the  colors.  An  example  of  this  effect  is  shown 
in  Figure  15  where  photographs  of  the  fringe  field  with  a  500A  wide 
filter  and  a  100A  wide  filter  are  displayed.  Very  good  fringes  are 
obtained  using  a  100A  wide  filter  centered  on  the  467  lA  xenon  line. 

4.  2  THERMAL  GRADIENTS  ABOUT  A  HORIZONTAL  CYLINDER 
IN  FREE  CONVECTION 

The  fringe  field  for  a  long  ceramic  cylinder  held  at  approxi¬ 
mately  900°C  in  air  with  free  convection  is  shown  in  Figure  16.  The 
horizontal  gradients  are  obtained  aligning  the  interferometer  so  that 
the  fringes  are  vertical.  Similarly  realigning  the  interferometer  the 
vertical  gradients  may  be  obtained  as  shown  in  the  upper  photograph. 
The  vertical  fringe  along  the  stagnation  line  shows  that  the  horizontal 
gradient  is  zero.  The  vertical  gradients  are  a  maximum  along  the 
same  line.  The  thermal  conductivity  could  be  determined  by  a  similar 
experiment  if  the  convective  flow  velocity  could  be  determined.  This 
experiment  would  avoid  many  of  the  radiative  and  convective  correc¬ 
tions  which  become  important  at  high  temperatures  and  pressures  in 
conventional  thermal  conductivity  measurements. 

4.  3  THEORY  OF  THERMAL  CONDUCTIVITY  MEASUREMENTS 
IN  THE  END  WALL  BOUNDARY  LAYER 

The  detailed  theory  of  the  Wollaston  prism  differential  inter¬ 
ferometer  is  given  in  References  6  and  26.  The  basic  theory  for  the 
thermal  conductivity  measurement  in  the  end  wall  boundary  layer  may 
be  understood  by  considering  the  two  rays  created  by  the  Wollaston 
prism  as  they  pass  through  the  shock  tube.  The  ray  traveling  through 
the  cooler  gas  travels  at  a  slightly  slower  speed  because  the  density 
of  the  gas  increases  as  it  cools.  The  two  rays  therefore  develop  a 
phase  difference  in  the  test  gas  given  by 


A  <f> 

2  TT 


± 

X 


(n+d-n-d> 


(20) 


where  £  is  the  length  of  the  test  section,  X  the  wavelength  of  the  light, 

n  .  and  n  ,  are  the  indexes  of  refraction  a  distance  y  +  —  from  the 
+d  -d  -  2 

end  wall.  The  theory  of  the  optics  gives  the  fringe  shift  As  in  terms 
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of  the  phase  shift  and  the  undisturbed  fringe  separation  s  by  the  relation 


As  =  s  =  s 

2tr 


n-d> 


or  finally  multiplying  and  dividing  by  the  ray  separation,  d,  we  obtain  a 
relation  between  the  density  (or  index  of  refraction)  gradient  and  the 
fringe  shift 


dp  _  _1  dn  _  _1  As  X.  _1 

dy  o  dy  a  s  £  d 


(22) 


where  a  is  the  Gladstone-Dale  constant  for  the  gas  in  question,  p  the 
density  of  the  gas  and  y  the  distance  from  the  end  wall. 

Equation  (22)  relates  the  measured  fringe  shift  As  to  the  density 
gradient.  The  relations  needed  to  obtain  thermal  conductivity  come 
from  the  boundary  layer  equations.  A  suitable  way  of  writing  the 
energy  equation  was  derived  in  last  year'  s  final  report  (Ref.  7).  The 
equation  is 

(y-yi)2  +  ^-rf  (y-yi)  +  c  (23) 


where  the  subscript  1  refers  to  the  reference  point  close  to  the  end  wall, 
Cpe  is  the  specific  heat,  and  X.e  the  thermal  conductivity.  The  similarity 
parameter  n  j  is  defined  by 


and  t  is  the  time  measured  from  the  time  that  the  shock  wave  is  reflected 
from  the  end  wall. 
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4.  4  APPLICATION  OF  DIFFERENTIAL  INTERFEROMETER 
TO  THE  WALL  BOUNDARY  LAYER 

The  steps  involved  in  applying  Equations  (22)  and  (23)  to  the  mea¬ 
surement  of  thermal  conductivity  are  (1)  derive  the  density  gradient  from 
the  measured  fringe  shift  on  the  interferogram  of  the  end  wall  boundary 
layer;  (2)  fit  the  density  gradient  to  the  quadratic  function  of  distance 
from  the  reference  point  (y^-y)^.  The  thermal  conductivity  is  given  by 
the  coefficient  of  the  quadratic  term  in  Eq.  (23).  Additionally,  the 
similarity  parameter  n  ^  may  be  estimated  to  within  a  few  percent  and 
the  linear  term  in  (23)  may  then  be  used  to  determine  thermal  conduc¬ 
tivity. 


There  are  several  important  points  about  Eq.  (23)  which  minimize 
experimental  errors  in  the  thermal  conductivity  measurements.  First, 
since  the  ratio  of  the  density  gradients  at  two  points  is  the  important 
property  in  Eq.  (23),  the  proportionality  factors  in  Eq.  (22)  do  not  need 
to  be  known,  that  is, 


\dy; 


In  addition,  the  magnification  in  the  direction  parallel  to  the  end  wall 
drops  out.  The  location  of  the  end  wall  is  also  not  necessary  because 
only  the  differences  in  distance  from  the  wall  occur  in  the  final  expres¬ 
sion  for  thermal  conductivity  (Eq.  23).  The  magnification  in  this  direc¬ 
tion  is  important  and  is  easily  determined  by  inserting  an  accurately 
known  cylinder  in  the  test  section  and  photographing  it.  Finally,  the 
accuracy  of  the  linearization  which  leads  to  Eq.  (23)  is  automatically 
determined  by  the  size  of  the  constant  (C)  which  approaches  zero  as 
the  deviation  from  free  stream  conditions  decreases  over  the  portion 
of  the  boundary  layer  in  the  schlieren  field. 

Some  measurements  of  the  fringe  displacement  on  a  typical  end 
wall  schlieren  interferogram  are  shown  in  Figure  17.  The  position  of 
the  end  wall  is  indicated  by  the  vertical  arrow.  The  data  are  measured 
points  on  the  fringe.  Taking  the  factor  of  two  magnification  on  the 
photograph  into  account  the  fringes  are  discernable  to  within  0.  1  mm 
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of  the  end  wall.  The  separation  of  ordinary  and  extraordinary  rays 
in  the  test  section  (the  rays  which  interfere  to  give  the  fringes)  is 
about  0.  047  mm,  much  smaller  than  the  boundary  layer.  The  total 
boundary  layer  is  3  mm  thick.  The  fringes  are  straight  lines  with 
~20%  slope  in  the  uniform  part  of  the  reflected  shock  region.  Closer 
to  the  end  wall  in  the  boundary  layer  the  fringes  bend  upward.  The 
total  temperature  change  in  the  part  of  the  boundary  layer  affecting 
the  fringes  is  about  3000  K.  This  means  that  half  of  the  temperature 
drop  is' in  the  outer  3  mm  of  boundary  layer  and  half  of  the  drop  is 
in  the  remaining  0.  1  mm.  The  temperature  gradient  near  the  end 
wall  is  very  steep,  while  the  outer  regions  show  small  temperature 
gradients.  This  is  the  basis  of  the  linearization  of  the  energy  equa¬ 
tion  to  give  the  useful  relation  Eq.  (23). 
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SECTION  V 

THERMAL,  CONDUCTIVITY  MEASUREMENTS  WITH  SCHLIEREN 

IN  TERFEROME  TER 


The  thermal  conductivity  of  nitrogen,  0.  25  argon- 0.  7^5  helium 
mixture,  and  air,  was  measured  at  temperatures  up  to  5000  K  and 
pressures  up  to  6  atm  using  the  schlieren  interferometer  on  the  end 
wall  of  the  shock  tube.  The  free  stream  temperature  and  pressure 
of  the  test  gas  behind  the  reflected  shock  were  determined  from 
incident  shock  speed  measurements  (Ref.  38).  The  transit  time 
determined  from  heat  transfer  gauges  106,  05  cm  and  14.  61  cm  from 
the  end  wall  was  used  to  determine  the  wave  speed.  In  addition,  an 
end  wall  heat  transfer  gauge  was  used. 

A  heat  transfer  trace  is  shown  in  Figure  17  along  with  a 
schlieren  interferometer  photograph  of  the  end  wall  boundary  layer. 

The  upper  heat  transfer  trace  is  the  side  wall  gauge  response,  the 
lower  trace  is  the  end  wall  gauge  response.  The  break  in  the  end 
wall  gauge  response  indicates  the  end  of  the  uniform  conditions  in 
the  gas  close  to  the  end  wall.  A  burst  of  rf  from  the  flash  tube  in¬ 
terrupting  uniform  end  wall  trace  indicates  when  the  schlieren  photo¬ 
graph  was  taken.  The  time  between  the  arrival  of  the  shock  wave  at 
the  end  wall  (the  step  in  the  end  wall  trace)  and  the  time  of  the  schlieren 
photograph  is  given  by  the  position  of  this  rf  pickup.  This  is  the  time 
which  should  be  used  in  Eq.  (23). 

Some  measurements  of  fringe  displacements  versus  distance 
from  the  end  wall  for  a  measurement  in  air  are  shown  in  Figure  18. 

The  ratio  of  -  is  averaged  over  a  number  of  fringes.  The  result- 

As  8 

ing  -  data  are  fit  to  a  quadratic  function  of  distance  as  described 

above  in  the  discussion  of  Eq.  (23).  The  resulting  thermal  conduc¬ 
tivity  values  are  tabulated  in  Table  4. 
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Table  4.  Thermal  conductivity  of  gases  measured  by  the 

schlieren  technique 


Temperature 

\ 

Nitrogen 

Pressure 

atm 

Thermal 
Conductivity 
mw/cm  °K 

2620 

2.  29 

2.  10 

4373 

4.  86 

3.  20 

4707 

5.44 

3.41 

5141 

6.  34 

4.  10 

Air 

2246 

1.  87 

1.87 

2258 

2.  04 

2.26 

2398 

2.  08 

2.  15 

2550 

1.  15 

2.  30 

2572 

0.  29 

2.25 

2999 

0.  30 

3.  00 

3157 

0.  57 

2.88 

4539 

5.  53 

3.  10 

0.  75  He  - 0.  25  Ar 

5571 

1.  24 

5.43 
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SECTION  VI 
SUMMARY 


Sound  absorption  and  sound  speed  measurements  were  shown  to 
be  a  useful  technique  for  measuring  gas  thermodynamic  and  transport 
properties  at  high  pressure.  The  sum  of  viscosity  and  thermal  con¬ 
ductivity  can  be  determined  from  sound  absorption  measurements  in  a 
wide  variety  of  pure  gases  and  gas  mixtures.  Measurements  of  sound 
absorption  in  argon- nitrogen  and  argon-helium  mixtures  were  made  at 
300°K  and  pressures  up  to  100  atm.  Good  agreement  between  theory 
and  experiment  was  obtained.  These  measurements  demonstrated  that 
losses  due  to  relaxation  and  diffusion  are  adequately  described  by  well- 
known  theory.  Therefore,  relaxation  and  diffusion  losses  may  be  de¬ 
termined  and  subtracted  from  the  total  absorption  to  obtain  the  sum  of 
viscosity  and  thermal  conductivity. 

The  sound  absorption  is  more  sensitive  to  the  density  dependence 
of  transport  properties  than  other  transport  property  measurements  be¬ 
cause  of  the  translational  (intrinsic  bulk  viscosity.  A  theoretical  relation 
between  bulk  viscosity  and  viscosity  would  be  useful.  Such  a  relationship 
can  probably  be  established  rigorously  in  the  near  future.  For  this  rea¬ 
son  ultrasonic  absorption  provides  support  for  long  range  theoretical 
developments  in  high  density  gas  transport  theory. 

Thermodynamic  properties  are  also  obtained  from  sound  speed 
measurements.  Thermodynamic  properties  of  mixtures  at  high  pressure 
are  not  generally  available;  therefore,  this  information  is  valuable  in  its 
own  right. 

Sound  speed  measurements  were  made  (with  accuracies  better 
than  0.  1%)  in  the  argon- nitrogen  and  argon-helium  systems  at  300  K  at 
pressures  up  to  400  atm.  Specific  heat  ratios  were  deduced  for  the 
argon-nitrogen  system.  Measurements  of  sound  speed  at  several  tem¬ 
peratures  and  at  pressures  above  100  atm  will  allow  a  complete  thermo¬ 
dynamic  description  of  these  gases  to  be  obtained. 

The  schlieren  interferometer  was  analyzed  for  thermal  conductivity 
measurements  at  temperatures  of  300  K  and  above.  Density  gradient 
profiles  were  measured  in  the  end  wall  boundary  layer  of  the  shock  tube. 
The  theory  of  end  wall  boundary  layer  was  used  to  derive  an  expression 
for  thermal  conductivity  in  terms  of  the  measured  density  gradient  pro¬ 
files.  Measurements  of  thermal  conductivity  of  nitrogen  and  air  were 
made  using  the  schlieren  interferometric  technique.  This  technique  is, 
in  the  authors'  opinion,  the  best  way  to  measure  thermal  conductivity 
at  temperatures  between  3000  K  and  8000  K. 
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FIG.  I:  PRESSURE-  TEMPERATURE  EFFECTS  ON  TRANSPORT 
PROPERTIES 
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FIG.  2:  CONTRIBUTION  OF  EXCESS  VISCOSITY  AND  BULK  VISCOSITY 
TERMS  TO  SOUND  ABSORPTION  IN  ARGON 
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Fig.  4.  Through  transmission  multiple  echo  pulses  in  air 
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FIG.  5 
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FIG.  6:  ABSORPTION  VS.  DENSITY  ARGON 


AE  DC-TR-69-78 


NORMALIZED  ABSORPTION  COEFFICIENT  -pax  I03  gm/cm 


PRESSURE  ATM 


FIG.  7s  ABSORPTION  VS.  DENSITY  IN  NITROGEN 
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FIO.  0!  ABSORPTION  VS.  DENSITY  50%  ARGON -50%  NITROGEN 
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FIO.  10:  ABSORPTION  VS.  DENSITY  -  52.6  ARGON  -  47.4  HELIUM 


AEDC-TR. 69-78 


8  ICE  I 
!_BATHJ 


FIG.  II:  ELECTRONICS  FOR  HIGH  PRESSURE  SOUND  ABSORPTION  AND  VELOCITY 
MEASUREMENTS 
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Upper  Trace  1  V/cm 
Lower  Trace  0.  5  V/cm 
2  |jl  sec/cm 
Argon  81.7  atm 


Fig.  12.  The  first  pulse  and  first  echo  used  for  absorption  and  velocity 
measurements. 
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FIG.  14:  SCHLIEREN  INTERFEROMETER 


100  A  filter 


500  A  filter 


Fig.  15.  Effect  of  monochromaticity  on  the  schlieren  interferometer  fringes 
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Vertical  gradients  T  =  900°C 
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Horizontal  gradients  T 
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Fig.  16.  Temperature  field  of  a  horizontal  cylinder  in  free  convection,  (x  5) 
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Fringe  Field  (x  2) 


Heat  Transfer  Trace 
200  p  sec/cm 


Fig.  17.  End  wall  boundary  layer  fringes  and  heat  transfer  gauge  response 
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FIG.  18:  FRINGE  FIELD  OF  END  WALL  BOUNDARY  LAYER 
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